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ABSTRACT 

We present a comparison between two observational and three theoretical 
mass functions for eight cosmological models suggested by the data from the 
recently completed BOOMERANG-98 and MAXIMA-1 cosmic microwave 
background (CMB) anisotropy experiments as well as peculiar velocities (PVs) 
and type la supernovae (SN) observations. The cosmological models have been 
proposed as the best fit models by several groups. We show that no model is 
in agreement with the abundances of X-ray clusters at ~ 1O 14 7 /i _1 M . On the 
other hand, we find that the BOOM+MAX+COBE:I, Refined Concordance and 
AMDM are in a good agreement with the abundances of optical clusters. The 
Pll and especially Concordance models predict a slightly lower abundances than 
observed at ~ W^h~ l M Q . The BOOM+MAX+COBE:II and PV+CMB+SN 
models predict a slightly higher abundances than observed at ~ 10 14 ' 9 /i _1 M Q . 
The nonflat MAXIMA-1 is in a fatal conflict with the observational cluster 
abundances and can be safely ruled out. 

Subject headings: cosmology:theory — cosmology:observation — galaxies: 
clusters: general — large-scale structure of universe 
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1. Introduction 

Recently, certain cosmo logical models have received a fairly strong observational boost. 
Several groups have used the new cosmic microwave background (CMB) data from the 
BOOMERANG-98 (the 1998 Balloon Observations Of Millimetric Extragalatic Radiation 
ANd Geophysics; |de Bernardis et al. 20001) and MAXIMA- 1 (The first overnight flight of 
the Millimeter Anisotropy experiment IMaging Array; [Hanany et al. 2000| ) anisotropy 
experiments to constrain cosmological parameters. Other groups combined the constraints 
from CMB with cosmological nucleosynthesis data, peculiar velocities (PVs) and type la 
supernovae (SN) observations. The values of cosmological parameters vary from one set to 
the next, but all of these models are in reasonable agreement with a flat Cold Dark Model 
(CDM) universe (f2 + I^a = 1) dominated by the vacuum energy except MAXIMA- 1 with 
matter density, Qq = 0.68 and vacuum energy density, Qa = 0.23 (Balbi et al. 2000). 

In this Letter, we compare the abundances of clusters of galaxies predicted by some 
popular cosmological models with observed abundances. The abundance of clusters has 
been shown to be one of the simplest but most effective cosmological tools for constraining 
the models of structure formation. It can place strong constraints on the parameters of 
cosmological models ( [Kaiser 1986 , Peebles, Daly, & Juskiewicz 1989, Simakov fc Shandarin 
1989 ), including the mass density in the universe (Qq) and the amplitude of the mass 



density fluctuations (erg) or, equivalently, the bias factor (b = l/a s ; |Evrard 1989| |Frenk et 



|al. 1990| ; |Henry fc Arnaud 1991| ; [Bahcall fc Cen 1992| ; |Lilje 1992| ; |Oukbir fc Blanchard 1992 



Kofman, Gnedin, & Bahcall 1993; White, Efsthatiou, & Frenk 1993; |Bond fc Myers 1996j 
Eke, Cole, & Frenk 1996; Mo, Jing, & White 1996; |Viana fc Liddle 1996| ; [Borgani et aT 



T997| ; |Henry 1997| ; |Pen 1998| ; |Postman 19991 ; |Verde et al. "2001] ; Pierpaoli, Scott, & White 



2001). 



The abundance of clusters and their evolution are quantified by the mass distribution 



-4- 



function. The theoretical derivation of the mass function of gravitationally bound objects 
has been pioneered by Press & Schechter (1974, hereafter PS). Despite various modifications 
that have been suggested recently (Cavaliere, Colafrancesco, & Scaramella 1991; Blanchard, 
Valls-Gabaud, & Mamon 1992; Monaco 1997(a,b); Audit, Teysser, & Alimi 1997; Lee & 
Shandarin 1998, hereafter LS; Sheth, Mo & Tormen 1999, hereafter SMT), it remains a 
viable model of the mass function and is widely used. 

In this Letter, we make use of three theoretical models suggested for the cosmological 
mass function: (i) the original PS mass function rips assuming the spherically symmetric 
collapse, (ii) the mass function n\ 3 that incorporates the anisotropic collapse as it is 
described by the Zel'dovich approximation (LS) and (iii) the mass function ust suggested 
by Sheth & Tormen (1999, hereafter ST) and later derived by SMT that takes into account 
both the anisotropic collapse and some nonlocal effects. Recently Jenkins et al. (2001) 
suggested fits to mass functions obtained in the "Hubble Volume" N-body simulations of 
some cosmological models. We have checked that using the fits by Jenkins et al. (2001) 
does not change the conclusions of this Letter. For comparison with observations we use 
mass functions obtained for cluster virial masses by Girardi et al. (1998) and Reiprich, 
Bohringer, & Schuecker (2000) 

Here we report the results for eight cosmological models. Among these, seven have 
recently been claimed as the best-fit models satisfying the data from CMB anisotropy 
experiments ( COBE Differential Microwave Radiometer, BOOMERANG-98 and MAXIMA- 
1) as well as from nucleosynthesis, large-scale structure and type la SN observations. 
These models are labeled Pll, BOOM+MAX+COBE:I, BOOM+MAX+ COBEJl, 
PV+CMB+SN, Refined Concordance, MAXIMA-1, and AMDM. We have included the 
Concordance model as a reference model since it is often referred to as the standard ACDM 
model. 
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None of the proponents of the best-fit models in our list have mentioned the explicit 
cluster abundance test. Rather some of them claimed that their models satisfy one of the 
many erg — relations reported in the literature, others even did not apply this test at all. 
We have noticed more than a dozen predictions of the as — relation in the literature, 
some of which are in conflict with the others. We believe our approach here to present the 
result of the cluster abundance test is more explicit. 

This Letter is organized as follows: in § 2 we briefly summarize the theoretical models 
of the cosmological mass functions, in § 3 we briefly describe the observational mass 
functions, in § 4 we outline the cosmological models, and, finally, in § 5 we report and 
discuss the results. 



2. Theoretical mass functions 

The cumulative mass function (cmf ) is the comoving number density of gravitationally 
bound objects of mass greater than M: N(> M) = n(M')dM' , where n(M)dM is the 
mass function of the collapsed objects with masses between M and M + dM. The PS model 
based on spherical collapse of overdense region in a smooth background predicts 

v 2 

np S (M) = F(p, a M ) v exp(-— ), (1) 

where v = 5 c /a M , F(p,a M ) = (2/n) 1/2 (p/M 2 ) ■ 

\d In &M/d In M\, and p is the mean matter density. The canonical value 5 C = 1.686 
corresponds to the spherical top-hat model in the Qo — 1 universe. Later it was shown 
that 5 C only weakly depends on the background cosmology ([Eke et al. 1996|) , and therefore 



we ignore it here. The rms density fluctuation (o"m) at the mass scale M is determined 
by the linear power spectrum a 2 M = l/2n 2 J °° dk k 2 P(k) W^ H (kR), where W T H{kR) 
is the Fourier transform of the top-hat window function. The mass M is related to R as 
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M = Air/3R 3 p. Although theoretically the most consistent approach requires the sharp 
k-space window function (see, e.g., |Bond et al. 1991|) , we use the top-hat window because 
it results in better fit to N-body simulations (see, e.g., ST). 

The A 3 -model suggested by LS is based on the non-spherical collapse as described by 
the Zel'dovich approximation. It assumes that a fluid particle belongs to gravitationally 
bound object after it experiences collapse along all three principle axes. In practice it 
has been approximated by imposing the condition A3 > A c on the smallest eigen value 
(A3 < A2 < Ai) calculated for the initial density field smoothed with the sharp fc-space 
filter corresponding to mass M. Comparisons with N-body simulations have shown that the 
threshold is A c = 0.37 (Lee & Shandarin 1999). The mass function in this model is given as 
(assuming A' = A c /<jm) 



n A3 (M) = ^Lf(p,<t m )\> 

24v27T 



20A'exp 



9A 



a 



2 



2tt(20A /2 - l) exp ( - ^)erfc( v / 2A' 
+ 3V3vrexp(^ — jerfc(^-j. (2) 

ST suggested a correction to the PS mass function resulting in better fit to N-body 
simulations (for a discussion of motivations see SMT) 

- ) v exp( 

a 



n ST (M) = F(p, a M ) A v 



1 + 1- 



av 

v exp( — — ). (3) 



The parameters A = 0.322, a = 0.707 and q = 0.3, chosen by ST, have been determined 
empirically from N-body simulation. At A = 1/2, a = 1.0 and q = 0, one finds ust = nps- 

The cosmological parameters enter the cosmological mass function via the shape and 
normalization of the linear power spectrum. One of the most accurate approximation of 
power spectrum fitting formula incorporating baryon density wass developed by Eisenstein 
& Hu (1998). Their formula has accuracy better then 5% for baryon fraction Qb/^o l ess 
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then 30%. The cosmological models discussed here predict baryon fraction less then 20%, 
therefore we have used the Eisenstein & Hu fits for the power spectrum. 

3. Observational mass functions 

The predictions of the theoretical models have been tested against the measurements 
of the virial mass functions in the N-body simulations (see, e.g., ST and references therein). 
Therefore, the theoretical mass functions must be compared with the observational virial 
mass functions. 

Girardi et al. (1998) provided the cumulative mass functions estimating the virial 
masses of clusters of richness R > — 1 and R > 1. Both practically coincide for 
M > 10 14M h~ 1 M Q (see Fig. 2 in Girardi et al. 1998). This mass function is shown by filled 
circles in Fig. 1 and 2. Reiprich et al. (2000) determined the cmf using X-ray flux-limited 
sample from ROSAT All-Sky Survey. They determined the masses from measured gas 
temperatures based on ASCA observations. In this Letter we use the mass function 
corresponding to r2oo which is usually referred to as the virial radius (open squares in Fig. 
1 and 2). At M < 1O 14 ' 8 /i _1 M the Girardi et al. mass function is significantly higher than 
that of Reiprich et al. 

It should be mentioned that the estimation of the masses is not a simple problem. For 
further discussion, see, e.g., Girardi et al.(1998), Reiprich et al.(2000), Pierpaoli et al.(2001), 
and references therein. In addition, there is no one-to-one correspondence between optically 
and X-ray-selected clusters. There are clusters found in both optical and X-ray surveys, but 
some optical clusters do not have counterparts in X-ray surveys and vice versa. There are 
some evidences suggesting that the fraction of X-ray clusters in a sample of optical clusters 
is smaller than the fraction of optical clusters in a sample of X-ray clusters. If confirmed by 
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further studies this means that some of optical clusters failed to become X-ray sources by 
some unknown reasons. However, this observation has been made for the ROSAT Optical 
X-ray Survey and must be taken with a great caution; it cannot be directly applied to any 
other surveys (M. Donahue 2001, private communication). Here we take both observational 
mass functions as they have been proposed by the authors without trying to resolve the 
discrepancies between them. 



4. Cosmological models 

In this Letter we discuss mostly flat cosmological models that are strongly motivated 
by the inflationary model of the universe (see e.g. Ostriker & Steinhardt 1995 and references 
therein). As an illustration we have included one open model (MAXIMA- 1 with Q tot = 0.91) 
advocated by Balbi et al. (2000) and one closed model (AMDM with Q tot = 1.06) advocated 
by Durrer & Novosyadlyj (2001). Although other groups (Valdarnini, Kahniashvili, & 
Novosyadlyj 1998 and |Primack fc Gross 200 ip have discussed AMDM type models, we have 



chosen only the above mentioned one for our comparison. The cosmological parameters 
have been obtained from observational data through likelihood analysis with various 
prior assumptions. These parameters (fib, ^cdm, &a, n s, h, a 8 ) from different models are 
presented in Table 1. In our notation, f2o = + Q c dm, spectral index n = n s + n t . In 
this letter, we have taken zero gravity wave contribution i.e. n t = with zero reionization. 
Among these models Pll, BOOM+MAX+ COBEI, BOOM+MAX+COBE:II, Concordance 
and MAXIMA-1 are COBE-normalized following the prescription of Bunn & White (1997). 
For other models we have followed the normalization suggested by the authors, models 



- 9- 



5. Summary 

We have compared the theoretical predictions of cluster abundance by several 
cosmological models with the observational mass functions determined by Girardi et 
al.(1998) (filled circles in Fig. 1,2) and Reiprich et al. (2000) (open squares in Fig. 1,2). 
In this Letter we make use of three theoretical mass functions rips, n\ s and nsr- It is 
worth stressing that in the range of masses (4 x 1O 14 /i _1 M < M < 3 x 10 15 /i _1 M Q ), the 
theoretical models differ one from another roughly less or similar to the error bars of both 
observational mass functions. 

At M < lO 14 8 /i _1 M no model can be reconciled with the Reiprich et al. (2000) X-ray 
mass function. On the other hand almost all models are in much better agreement with 
the Girardi et al. (1998) optical mass function. Thus, the resolution/explanation of the 
discrepancies between optical and X-ray mass functions becomes crucial for the well being 
of all models in question. 

As far as the optical mass function is concerned the Refined Concordance, 
BOOM+MAX+ COBE.l, and AMDM models show a reasonable agreement with 
observations. The Pll and especially Concordance models predict a slightly 
lower abundances than observed at ~ 10 14 ' 6 /i _1 M Q . On the other hand, the 
BOOM+MAX+COBE:II and PV+CMB+SN models predict a slightly higher abundances 
than observed at ~ 10 w - 9 h~ 1 M ( 7 ) . The MAXIMA-1 model seems to be safely ruled out by 
the data on cluster abundances. 

A similar comparison using the sharp k-space filter for evaluation of o"m, which is better 
justified for the PS mass function (|Bond et al. 1991) , showed that all three theoretical mass 



function are systematically higher than that for the top-hat filter. The sharp k-space filter 
approach improves the agreement with observations for the Pll and Concordance models 
and makes it worse for the BOOM+MAX+COBE:II and PV+CMB+SN models. The 
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conclusions for other models did not change much. 
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Table 1 

Parameters of the Cosmological Models 



Models 


Parameters 


Reference 




^cdm 




n s 


h 


^8 


Pll 


0.045 


0.255 


0.7 


0.95 


0.82 


0.92 


Lange et al.2001 


BOOM+MAX+CttBE:I 


0.045 


0.255 


0.7 


0.975 


0.82 


0.97 


Jaffe et al.2000 


BOOM+MAX+ COBEJl 


0.036 


0.314 


0.65 


0.95 


0.80 


1.06 


Hu et al.2001 


PV+CMB+SN 


0.035 


0.245 


0.72 


1.0 


0.74 


1.17 


Bridle et al.2001 


Concordance 


0.03 


0.27 


0.7 


1.0 


0.68 


0.85 


Ostriker & Steinhardtl995 


Refined Concordance 


0.05 


0.33 


0.62 


0.91 


0.63 


0.83 


Tegmark et al.2001 


MAXIMA-1 (n tot = 0.91) 


0.07 


0.61 


0.23 


1.0 


0.60 


1.05 


Balbi et al.2000 


AMDM (Q tot = 1.06) 


0.037 


0.303 


0.69 


1.02 


0.71 


0.92 


Durrer & Novosyadlyj 2001 



This manuscript was prepared with the AAS IATj^X macros v4.0. 
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Figure Captions 

Fig. 1. Observational cmfs measured for virial mass are compared with different 
theoretical predictions: (a) Pll, (b) BOOM+MAX+ COBE: I, (c) B00M+MAX+COBE: 
II and (d) PV+CMB+SN model. The short dash line is nps, long dash line is n\ 3 and solid 
line is nsr] the filled circles are the observational data points corresponding to virial masses 
determined by Girardi et al. (1998) and the open squares are those determined by Reiprich 
et al (2000). The open triangle is the value of the cmf for masses estimated within the 
1.5/& -1 Mpc radius given by Girardi et al. The error bars are in la limit along the vertical 
direction. Horizontal bars indicate the bin size. 

Fig. 2. Same as Fig. 1 but with different models: (a) Concordance, (b) Refined 
Concordance, (c) MAXIMA-1 and (d) AMDM model. 



-16- 




Fig. 1.- 



-17- 




14.6 14.8 15 15.2 15.4 14.6 14.8 15 15.2 15.4 

iog 10 [ M/(h-» M e ) ] 



Fig. 2.- 



